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PHOTON BUBBLES IN YOUNG MASSIVE STARS 
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ABSTRACT 

Spectroscopic studies indicate that gas in the photospheres of young O stars moves at speeds up to the sound 
speed. We show, using two-dimensional radiation MHD calculations and results from a local linear analysis, 
that the motions may be due to photon bubble instability if young O stars have magnetic fields. 
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1. INTRODUCTION 

Young high-mass stars may trigger or terminate 
nearby star and planet formation through ionizing 
radiation and winds launched from their surfaces. 
The properties of the radiation and winds depend 
on the structure of the stellar surface layers. Stars of 
more than 15 Solar masses reac h the main sequence 
while still accreting material ()Yorke Sonnhalterl 
I2OO2I) . High- mass main sequence stars are thought 
to have stably stratified outer layers, yet gas motions 
are present in the atmospheres of these stars. 

Absorption lines formed in the photospheres of O 
stars are broader than expected based on the tem- 
perature, pressure and stellar rotation. If the rota- 
tion axes are randomly oriented with respect to our 
line of sight, some stars will be viewed pole-on and 
show no rotational broadening. However samples 
of O stars appe ar to contain no such pole-on cases 
lS1etteba,]Jh956tlConti Ebbet jlToTTi IPenn ^11993 
iHowarth et al.lll997l) . The minimum apparent equa- 
torial rotation speed is a function of spectral type 
and luminosity. Among main-sequence O stars in the 
Small Magellanic Cloud cluster NGC 346, the excess 
broadening of UV metal lines ranges from 25 km s~^ 
at 02, similar to the speed of sound at the ph oto- 
sphere, to 5 km s"^ at 09.5 (jSouret et al.ll2003l) . 

Processes that might be responsible for the addi- 
tional line broadening include 

1. Vertical velocity g radients in the ste llar wind ac- 
celeration region ()Kudritzkilll992l) : Minimum 
line widths are similar in stars with and without 
strong outfiows. Also, stellar atmosphere plus 
wind calculations including a variation in veloc- 
ity through the acceleration region produce pho- 
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tospheric lin es similar to those in hydrostatic 
calculations (|Bouret et al.ll2003l) . 

2. Global non-radial pulsations: Low levels of spec- 
tral li ne shape variat i ons in main- sequence O 
stars ([Pullerton et al.l Il996l) suggest that any 
global non-radial modes are either of low am- 
plitude or high order. An unresolved issue is 
how such modes might be excited. 

3. Strange modes: Instability occurs only for 
masses greater than 80 M© in zero-age main 
sequence stars with meta llicity Z = 0.02 
(|Glatzel Kiriakidislll993l) . while excess pho- 
tospheric line widths are observed in main- 
sequence stars with masses down to 20 M©. 

4. Shaviv modes: The instabilities are found only 
at lu minosities grea ter than half the Eddington 
value iShavivll200ll) . while excess line broaden- 
ing occurs in stars with one-tenth the Eddington 
luminosity. 

Here we explore the possibility that a fifth 
mechanism, ph oton bubble instability (lAronslll99i 
lGammijll998j) . may lead to small-scale motions in 
the surface layers of O stars. A general local WKB 
analysis of linear disturbances in optically-thick radi- 
ating atmospheres i ndicates instability und er a wide 
range of conditions (|Blaes Socratesi2003l) . Photon 
bubbles are driven by the radiative fiux (Figure^, 
and grow if (1) Rosseland mean optical depth per 
wavelength is greater than unity, (2) a magnetic field 
is present, and (3) the radiative fiux exceeds a criti- 
cal value that in radiation-supported atmospheres is 
approximately the gas sound speed times the sum of 
gas and photon energy densities. The surface layers 
of O stars generally satisfy the first and third criteria, 
but until recently were thought to be unmagnetized. 
However 0^ Ori C, the illuminating star of the Orion 



1 



2 



TURNER ET AL. 



Gravity 



Magnetic 
FieU 



4 



Radiation Flux 



Gravity 



Magnetic 
Field 



Radiation Flux 



Fig. 1. How photon bubbles work. Toy: A patch of at- 
mosphere (gray square) is initiahy in hydrostatic balance. 
When density is reduced slightly in a small elongated re- 
gion, photons diffuse more easily up the long axis. The 
perturbed radiation flux exerts an extra force with a com- 
ponent along the magnetic field. Gas is driven along the 
field to the right, out of the region of low density, and the 
perturbation grows. Bottom: As fresh gas enters the re- 
gion of low density from the left, the pattern propagates 
to the left at about the gas sound speed. 



Nebula and one of the nearest young O stars, shows 
polarimetric variations in photospheric lines with ro- 
tational phase, consistent with a 1.1 k G dipole field 
inclin ed 42° from the rotation axis ([Ponati et alJ 
Also, narrow X-ray emission lines indicate 
some gas near the star moves more slowly than the 
wind. The line ratios show that most of the X-ray 
flux is thermal emission from plasma hotter than 15 
million K, too hot to arise in shocks in the wind, 
and as hot as co ronae in sonie mag netically active 
lower-mass stars ([Schulz et al J 120031) . Overall, there 
is good evidence for a magnetic field in 0^ Ori C. 

Photon bubbles, unlike strange modes, can grow 
in scattering atmospheres. They differ from Shaviv 
instabilities in that they require magnetic fields and 
may grow when the luminosity is substantially less 
than the Eddington value. 



2. METHODS 

A small patch of the surface layers of the star is 
followed in Gartesian coordinates {x^y^z). The x- 
axis is horizontal, the z-axis vertical, and symmetry 
is assumed along the ^-axis. The domain is 10^^ cm 
on a side, or about 1.5% of the radius of 6^ Ori G, and 
divided into 32 x 32 zones. The top of the domain 
is placed at the photosphere. The side boundaries 
are periodic, and the upper and lower boundaries are 
transparent walls through which radiation may pass, 
but gas and magnetic fields do not. The tempera- 
ture of the lower boundary is fixed at its initial value, 
while the radiative flux is continuous across the up- 
per boundary. The equations of radiation MHD, 
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2OO3I) are solved u sing the ZEUS code 
([stone Normanlll992allb[) with its fiux -limited ra- 
diation diffusion module ([Turner Stone 2001.) . A 
uniform gravity g is included. Rotation is ignored, as 
photon bubbles grow to saturation within a typical 
stellar rotation period. Emission proportional to the 
blackbody rate B is assumed in equations 01 and ^ 
Temperatures are calculated using a mean mass 
per particle of 0.6 protons. The total flux-mean 
opacity xf includes contributions from electron 
scattering 0.34 cm^ g~^ and bound- free and free- free 
processes 2 x lO^^p^'^e"^'^ cm^ g~^, where e is the 
gas internal energy per unit volume and p is the 
mass density. The Planck-mean absorption opacity 
tvp due to the bound- free and free- free effects is 
7.4 X 10^^/-^e-^-^ cm^ g-^ The radiative flux 
F = — —WE is calculated using; a limiter A that 

XFP ^ 

equals ^ in optically-thick regions, and is reduced in 
regions of low optical depth so that the flux does not 
exceed the product of the radiation energy densit y 
E and light speed c ([Levermore PomraninJl98lh . 
The set of equations is closed using an ideal gas 
equation of state p = (7 — l)e with 7 = 5/3. 
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Initial conditions are constructed by integrating 
the equations of hydrostatic equihbrium and radia- 
tive flux conservation from the photosphere down 
into the interior. Three initial states are considered, 
spanning the range of effective temperatures found in 
O stars. The photospheric temperatures in Kelvins 
and densities in g cm~^ are (A) 50 000,1.8 x 10~^, 

(B) 40 000,2.2 X 10-^ and (C) 30 000,3.0 x 10"^ 
Surface gravity 10^"^ cm s~^ is used in all cases. The 
Eddington ratios between the last two terms in equa- 
tion[2lat the photospheres are (A) 0.65, (B) 0.28, and 

(C) 0.11. We study the growth of photon bubbles at 
the stellar magnetic equator by beginning with a uni- 
form, horizontal magnetic field of 1 kG. Radiation, 
gas and magnetic pressures at the photosphere are in 
ratios (A) 0.40 : 0.31 : 1, (B) 0.16 : 0.31 : 1, and (C) 
0.05 : 0.31 : 1. Random zone-to-zone initial density 
perturbations with amplitude 0.1% are added only 
in the middle half of the domain height, to reduce 
any effects of the boundaries on the linear growth. 

3. RESULTS 

Each calculation passes through three stages. 
First, the perturbations cause small damped verti- 
cal oscillations with period equal to the gas sound 
crossing time and amplitude 0.1 km s~^. Next, 
linear eigenmodes appear and grow exponentially. 
The fastest-growing modes are gas sound waves with 
wavefronts inclined 30 to 45 degrees from vertical, 
and group velocity roughly parallel to the magnetic 
field. Finally, the growth saturates and amplitudes 
vary little over many gas sound crossing times. 

3.1. Linear Growth 

The time histories of the horizontal kinetic en- 
ergies are shown in Figure [2l Linear growth rates 
increase with photospheric temperature. In the 
50 000 K case (A), the fastest mode has 3 wave- 
lengths in the domain width and 2 in the height, cor- 
responding to wavelength 2.8 x 10^ cm and wavevec- 
tor 34° from horizontal. The power in the mode, 
measured using the Fourier transform of the horizon- 
tal velocity pattern, grows e-fold every 3000 seconds. 
The growth rate is 2/3 that obtained by solving the 
linear dispersion relation, Blaes & Socrates (2003) 
equation 49, for conditions at the center of the com- 
putational domain. The slower growth in the simu- 
lation may be due to the modest grid resolution. 

3.2. Saturation 

The disturbances reach larger amplitudes in the 
cases with hotter photospheres. The largest horizon- 
tal velocities are 30 km s~^ in the 50 000 K case (A), 
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Fig. 2. Domain- averaged kinetic energies in horizon- 
tal motions versus time, in the calculations with photo- 
spheric temperatures 50 000 (top), 40 000 (middle), and 
30 000 K (bottom) . The upper horizontal line indicates 
the energy density in the initial magnetic field. 

12 km s~^ in the 40 000 K case (B), and 6 km s~^ in 
the 30 000 K case (C). In the hottest case (A), the 
speeds are similar to the isothermal gas sound speed 
of 26 km s~^ at the photosphere, and weak shocks 
are present. The structure resembles that proposed 
by Begelman (2001), but with crossed shocks due to 
the initial horizontal symmetry. Gas oscillates back 
and forth along magnetic field lines, experiencing re- 
peated compression and expansion. Densities vary 
about their initial values by up to a factor three, 
while radiation energy remains almost uniform due 
to rapid diffusion (Figure|3|). Radiation escapes more 
easily through the inhomogeneous atmosphere than 
through the original hydrostatic structure, and the 
time-averaged radiative flux is 10% greater than ini- 
tially. The density variations in the 40 000 K case 
(B) are 50%, and in the 30 000 K case (C) are 20%. 
The magnetic fields bend 30° up and down in the 
hottest case, and remain close to horizontal in the 
two cooler cases. 

4. CONCLUSIONS 

Photon bubble instability is present in the 
surface layers of magnetized main-sequence O stars 
according to the results of a local linear WKB 
analysis. We use two-dimensional radiation MHD 
calculations with horizontal magnetic fields and 
closed vertical boundaries to show that the insta- 
bility can lead to small-scale movements of gas in 
the surface layers, with velocities about equal to the 
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Fig. 3. Snapshot of saturated photon bubble insta- 
bihty in the calculation with photospheric temperature 
50 000 K, at 56 000 seconds. The initial photosphere 
Ues along the top boundary, and the domain height and 
width are 10^° cm. Top: Density on a gray scale loga- 
rithmic between 1.2x 10~^ (black) and 1.3x 10~^ g cm~^ 
(white). Velocities are shown by arrows, the longest cor- 
responding to 29 km s~^. Bottom: Radiation pressure on 
a gray scale linear between 6 600 and 210 000 dyn cm~^. 
Magnetic fields are shown by arrows, the longest 1500 G. 

observed linewidth excesses. In a calculation with 
parameters appropriate for an early O star, photon 
bubbles result in horizontal density variations near 
the photosphere. Radiation diffuses more rapidly 
through the inhomogeneities than through the 
hydrostatic atmosphere with the same column 
depth, and the radiative flux is 10% greater than 
that calculated assuming hydrostatic equilibrium. 
An enhanced flux may affect measurements of 
stellar parameters. Issues remaining to be resolved 



include the dependence on the local magnetic field 
orientation and the shapes of the spectral lines ex- 
pected from emission integrated over the stellar disk. 

Photon bubbles reaching amplitudes similar to 
those in the hottest case we considered may eject 
blobs of gas through the photosphere, and could lead 
to density variations in the winds of O stars. The 
motions might transfer energy into an initially weak 
magnetic field. Additional radiation MHD calcula- 
tions indicate the saturation amplitude generally in- 
creases with decreasing surface gravity, so photon 
bubbles may reach larger amplitudes in O giants 
and supergiants. Solutions of the dispersion rela- 
tion show that photon bubbles may be present also 
in accretion disks around young massive stars. 

N. J. T. was supported by a National Research 
Council fellowship. The work was carried out in part 
at the Jet Propulsion Laboratory, California Insti- 
tute of Technology. 
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